args <- commandArgs(TRUE)
geno<-read.table(args[1],sep="\t",head=T)
gnp<-read.table(args[2],sep="\t",head=T)
geno1<-as.vector(as.numeric(geno[1,]))
gnp1<-as.vector(as.numeric(gnp[1,]))
 gnp2<-as.vector(as.numeric(gnp[2,]))
gnp3<-as.vector(as.numeric( gnp[3,]))
d<- data.frame(gnp1,gnp2,gnp3,geno1)

glm1 <- glm(geno1 ~ gnp1+gnp2+gnp3, data=d,family=gaussian(link="log"))
glm1
pred <- predict(glm1, newdata=data.frame(gnp4), type="response")
pred

